The integrated use of electricity and natural gas has captured great attention over recent years, mainly due to the high efficiency and economic considerations. According to the energy hub design and operation, which allows using different energy carriers, it has turned into a critical topic. This paper develops a two-stage stochastic model for energy hub planning and operation. The uncertainties of the problem have arisen from the electric, heating, and cooling load demand forecasts, besides the intermittent output of the solar photovoltaic (PV) system. The scenarios of the uncertain parameters are generated using the Monte-Carlo simulation (MCS), and the backward scenario reduction technique is used to alleviate the number of generated scenarios. Furthermore, this paper investigates the effectiveness of demand response programs (DRPs). The presented model includes two stages, where at the first stage, the optimal energy hub design is carried out utilizing the particle swarm optimization (PSO) algorithm. In this respect, the capacity of the candidate assets has been considered continuous, enabling the planning entity to precisely design the energy hub. The problem of the optimal energy hub operation is introduced at the second stage of the model formulated as mixed-integer non-linear programming (MINLP). The proposed framework is simulated using a typical energy hub to verify its effectiveness and efficiency.